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Abstract. 

Starting from the three-dimensional (3D) time-dependent nonlinear Gross- 
Pitaevskii equation for a Bose-Einstein condensate (BEC) and density functional (DF) 
equation for a Fermi superfluid at the unitarity and Bardecn-Coopcr-Schrieffer (BCS) 
limits, we derive effective one- (ID) and two-dimensional (2D) mean-field equations, 
respectively, for the dynamics of a trapped cigar- and disk-shaped BEC and Fermi 
superfluid by using the adiabatic approximation. The reduced ID and 2D equations 
for a cigar- and disk-shaped Fermi superfluid have simple analytic nonlinear terms and 
at unitarity produce results for stationary properties and non-stationary breathing 
oscillation and free expansion in excellent agreement with the solution of the full 3D 
equation. 
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1. Introduction 

Many of the stationary and non-stationary properties of a trapped Bose-Einstein 
condensate (BEC) (or a Bose superfluid) can be satisfactorily explained by the mean- 
field Gross-Pitaevskii (GP) equation pQ. However, there is no such mean-field equation 
in configuration space for the Fermi superfluid even after the experimental realization 
of the crossover from the Bardeen-Cooper-Schrieffer (BCS) limit through unitarity to 
the Bose limit in a trapped Fermi superfluid [21 [31 HI [5] . 

To study the properties of stationary and non-stationary states of a trapped Fermi 
superfluid with equal number of paired spin-up and -down fermions, Adhikari and 
Salashnich developed a Galilei-invariant nonlinear density functional (DF) equation valid 
from the weak-coupling BCS limit to unitarity [61 [7] . The solution of the DF equation 
of [6] yielded results for energy of a trapped Fermi superfluid in close agreement with 
those obtained from Monte Carlo calculations [8] not only in the BCS and unitarity 
limits but also along the BCS-unitarity crossover 

In actual experiments, an axially symmetric [U [3j [9] , rather than a spherical, trap 
is usually employed for the confinement of the BEC or Fermi superfluid. In many 
situations of actual interest, the trap has an extreme geometry [9], e.g., very strong 
radial or axial confinement. Consequently, the superfluid is formed in the shape of a 
cigar or a disk and the solution of the nonlinear equation in such cases deserves special 
attention. The cigar-shaped BEC or Fermi superfluid is quasi one dimensional (ID) and 
the reduction of the full three-dimensional (3D) equation to accurate ID form would 
be of advantage. Similarly, it would be beneficial to reduce the full 3D equation to 2D 
form for the description of a disk-shaped BEC and Fermi superfluid. 

We propose a simple scheme for the 3D-1D and 3D-2D reduction of the GP equation 
for a BEC and the DF equation for a Fermi superfluid at the BCS and unitarity 
limits. The 3D- ID reduction for a cigar-shaped BEC or a Fermi superfluid is done 
in the adiabatic approximation which assumes the essential dynamics to be confined in 
the axial direction with the radial degrees of freedom adjusting instantaneously to the 
minimum energy equilibrium configuration compatible with the axial dynamics [TU| [TT]. 
Under this approximation the 3D wave function factorizes into an explicitly time- 
dependent axial and time-independent radial parts which allows for a formal reduction 
of the original time-dependent 3D nonlinear equation into a time-dependent ID and 
an auxiliary time-independent 2D equations. The same is also true for a disk-shaped 
BEC or a Fermi superfluid with the role of time-dependent ID and auxiliary time- 
independent 2D equations interchanged, e.g., the time-dependent equation is now 2D 
and the auxiliary equation ID in nature. 

First we illustrate the present scheme through an application to a BEC described 
by the GP equation, where there are already several schemes pUl HU C21 [13], [HI [151 021 
ITTl [T8| [T9] for 3D-1D and 3D-2D reductions. Among these, the reduction schemes of 
Salasnich et al. (SPR) [TO, E2j and Munoz Mateo et al. (MMD) [TT1[T3] are the simplest 
to implement and have been shown to be the most efficient [UJ . We compare the present 
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scheme with those of SPR [TQl H2] and MMD [HI E3] for a BEC. Recently, Adhikari and 
Salasnich (AS) [20] suggested one scheme for such 3D-1D and 3D-2D reductions of the 
3D DF equation for a Fermi superfluid at the BCS and unitarity limits. We apply our 
scheme to a Fermi superfluid and obtain distinct 3D- ID and 3D-2D reductions. 

In section [2] we illustrate the 3D- ID and 3D-2D reduction schemes in the case 
of the GP equation for a BEC. In section [3] we consider the reduction schemes for the 
nonlinear DF equation of a Fermi superfluid at unitarity and BCS limits. In this case we 
reduce the 3D DF equation to ID (for cigar-shaped superfluid) and 2D (for disk-shaped 
superfluid) forms with analytic nonlinear terms. We show that the reduced ID and 2D 
equations for cigar- and disk-shaped Fermi superfluids are very effective for describing 
the stationary states and non-stationary breathing oscillation and free expansion of the 
Fermi superfluid at unitarity in close agreement with the solution of the full 3D DF 
equation. Finally, in section 4 we give a summary and some concluding remarks. 



2. 3D-1D and 3D-2D Reductions of the GP equation 

We first apply our approach to a BEC described by the GP equation. The GP equation 
for N bosons of mass m is written as [1] 



-%h— - + V(r) H *W 

dt 2m K J m 1 v n 



tt(f,t) = 0, (1) 



with a the atomic scattering length, normalization / |\I/| 2 <i 3 ? = 1, V(f) = mu 2 (p 2 + 
X 2 z 2 )/2 the harmonic trapping potential with frequencies u and Xuj in radial (p) and 
axial (z) directions (f = p, z). Employing dimensionless harmonic oscillator units 
t = iu, r = r/a p , z = z/a p ,p = p/a p ,a = a/a p ,ip = ^a^ 2 ,a p = ^Jh/mu, ((TJ can be 
written as 

72 ^2 i \2^2 

V(r,t) = 0, (2) 



. d V r 2 p 2 + A V 

-i- + h AnaN\ 

dt 2 2 1 



with normalization / |?/;| 2 <i 3 r = 1. 

In the limits of very small (cigar-shaped trap) and very large (disk-shaped trap) 
A, the length scales in the axial and radial directions are very different. Consequently, 
the correlations between these two directions can be neglected and the condensate wave 
function could be factorized in variables p and z pUJ [H]. In the case of cigar-shaped 
traps the dynamics takes place in the axial direction. The opposite happens in case of 
a disk-shaped trap. For a stationary solution ^(r) of ([2]) in a cigar-shaped trap one can 
define a linear density <fi 2 {z) = / d 2 pif) 2 (r). Similarly, for a disk-shaped trap one can 
define a radial density v ?2 (p) = / dz%fj 2 {r). 

2.1. 3D-1D reduction for a cigar-shaped BEC 

For a cigar-shaped trap, the above consideration leads to the factorization [11] 

ip(r, t) = <p(p, m(z, t))<j)(z, t), (3) 
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where linear density n\ is defined as ni(z,t) = N\(p(z,t)\ 2 = Nfd 2 p\ip\ 2 and 
normalizations / d 2 p\ip(p, ni)\ 2 = f dz\<j)(z,t)\ 2 = 1. We assume that the function 
(p(p, ni(z, t)) has no explicit t or z dependence and hence these derivatives do not act on 
this function. This is a good approximation, in general, as we find in numerical studies. 
The substitution of flU in O) then leads to 



<p(p,ni 
Multiplying 



.d Id 2 
l di + 2lh 2 



-X 2 z 2 
2 



<f){z,t)= <f>(z,t) 



-V 

2 p 



—p 2 + 47rani|9?| 



ip(p,n x ). (4) 



by ip*(p,ni) and integrating in p this equation can be rewritten as [TT 
.d Id 2 1 



2^2 



dt 



2dz 2 2 

;2 



=0, 



1 

2 ' 9 ' 2' 

where p p {ni) is the chemical potential emerging from the 2D GP equation: 



-p 2 + 4irani\(p\ 



1 2 1 2 

2 p 2^ 



4:7iani\ip\ 2 — p p (ni] 



o. 



(5) 
(6) 

(7) 



We have decoupled the essential axial (z) and nonessential radial (p) degrees of freedom. 
The solution of the time-independent radial GP equation Q leads to the chemical 
potential p p {ni) given by ([6]), which is the nonlinear term of the axial GP equation (j^J). 

The form of the chemical potential p p (ni) of ([7]) is known in the small and large N 
limits. In the small N weak-coupling limit, the wave function can be approximated by 
the following normalized Gaussian form 



<P(P, rii) = exp[-p 2 / (2a 2 )]/ (vVir), 
where a is the width. With this wave form the chemical potential of 



becomes 



a 

T 



l 



2a 5 



1an\ 



a- 



(8) 



(9) 



In the large Thomas- Fermi (TF) limit, as an x — > oo, the kinetic energy gradient 
operator in ([7]) can be neglected and this equation has analytic solution. The 
normalization condition of the TF wave function leads to [TQl ITT] 

P P (ni) 



y/4an~i. 



(10) 



We suggest the following simple interpolation formula for p p (ni) valid for small to 
large an\ incorporating the limiting values (Q and fTTOl) 



p p {ni) 



1 



2a 2 



or 
T 



a 4 + 4ani 



to be used in dSJ), where a is taken as a fixed constant for all n\. In the weak-coupling 
ari\ — >• limit, (J7J) reduces to the Schrodinger equation for a 2D harmonic oscillator with 



the exact solution (Jo]) with a = 1. Then (TTTj) is a good approximation to (Jo]) for a = 1. 
For slightly larger values of an x fJTT]) continues to be a good approximation to (JBj), but 
with a slightly smaller value of a. Motivated by this, we take a slightly smaller value 
of a in ( TTTT) . For large ani, ( TTTT) has the correct TF limit independent of the value of a 
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employed. This flexibility in introducing a width a ~ 1 (slightly different from a = 1) in 
(Hip will be fundamental in making the ID model equation ([5]) a faithful approximation 
to the 3D GP equation ([2]) for a cigar-shaped condensate for all an\. 

By construction, approximation (TTTT) satisfies the weak-coupling and TF limits (Q 
and ( fTOl) . respectively, for small and large an\. The approximation of MMD is [H] 



//p(ni) = Vl + 4ani, (12) 





- 4ani, 






i + 


3arii 



^ (ni) = v /l + 2an 1 ' (13) 

Here we calculate the chemical potential p, p (ni) of the three approaches and 
compare with the precise results for fi p (ni) obtained from a numerical solution of (J7J). 
(All numerical results presented in this paper are obtained with the imaginary-time 
propagation scheme after a discretization by the Crank-Nicholson method using the 
FORTRAN programs provided in [22], the details of which are described there. The 
numerical simulations for the dynamical breathing oscillation and free expansion for 
a Fermi superfluid at unitarity reported in Sees. IIIC and HID, respectively, were 
performed with the real-time propagation scheme after a discretization by the Crank- 
Nicholson method.) Our finding is exhibited in tabledJfor different ari\ and for a = 0.985 
together with those obtained from the MMD [TT] and SPR [TU] schemes. 



Table 1. Chemical potential /z p (ni) of © for different an\ obtained from an accurate 
numerical solution of ©, and from (TTT]) (Present, a = 0.985), HU) (MMD) [TT], and 

PI (spr) pro]. 



an i 


Numerical 


MMD 


SPR 


Present 





1 


1 


1 


1.00045 


0.2 


1.346427 


1.34164 


1.35225 


1.34983 


1 


2.257135 


2.23607 


2.30940 


2.25314 


10 


6.432456 


6.40312 


6.76475 


6.42877 


100 


20.04320 


20.0250 


21.2309 


20.0538 



2.2. 3D-2D reduction for a disk-shaped BEC 

For a disk-shaped trap the adiabatic approximation leads to the factorization [101 [11] 

i[){r,t) = <p(p,t)(f>(z,n2(p,t)), (14) 
where the surface density n 2 is defined as n 2 (p,t) = N\<p(p, t)\ 2 = N f dz\ip\ 2 and 



normalizations / d 2 p\(p(p, t)\ 2 
leads to 

1 



■ 9 1^ 2 
*8t + 2 V » 



2 P 



J dz\<p(z,n 2 )\ 2 



1. The substitution of ffl~4l) in 



l^ + l^ + AnanM 2 



<t>(z,n 2 )(15) 
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Multiplying (11 5p by (f)*(z, n 2 ) and integrating in z, this equation can be rewritten as the 
set of equations 

1 n 



■ 9 1^9 



2 P 



<p(p,t) = 0, 



1 d 2 \ 2 z 2 
+ 



2dz 2 



2 

dzc 



Airan 2 1 



2 - Hz{n 2 ) 
\ 2 z 2 



,z,n 2 ) 



4iran 2 



(16) 
(17) 
(18) 



It is convenient to introduce scaled variables z 
a z = Jljx, when f fT7|) and ffTB"]) become 

1 d 2 



z/a z ,(f> = y/a z (f), and \i z = fi z a 2 with 



z 2 



[z,n 2/ 



Hz{n 2 ) 



dzcfi* 



1 d 2 z 2 

--— — H h 4iiaa z n 2 \(f)\ 2 

2dz 2 2 m 



(19) 
(20) 



The form of the chemical potential j2 z (n 2 ) of ffl9l is known in the small and large 
N limits. In the small N weak-coupling limit the wave function can be approximated 
by the following normalized Gaussian form [21] 

4>(z,n 2 ) = exp[-z 2 /(2^)]/(^) 1/4 , (21) 
where 77 is the width. With this wave form the chemical potential of ( 1201) becomes 

r + TT + 2aa zn 2 • 

4 Ar] z J rj 



(22) 



AM n 2j = 

In the large N limit the normalization condition of the TF wave function leads to [11] 

pzM = (3Tiaa z n 2 /V2) 2 / 3 . (23) 

For a disk-shaped BEC we suggest the following simple interpolation formula for 
^2(^2) incorporating the limiting values (122!) and (123]) 



4 



7T?7 



2\ 3/2 



4 



37raa 2 n2 



2/3 



(24) 



valid for all aa z n 2 , where 77 is taken as a fixed constant for all aa z n 2 . The flexibility in 
introducing a width 77 slightly different from 1 in (1241) will make the 2D model equation 
( TIB"]) an accurate approximation to the 3D GP equation ([2]) for a disk-shaped condensate 
for all aa 2 z n 2 . 

Instead of taking 77 as a constant, SPR [10] solve (fl9~l) variationally [10] [21] with the 
Gaussian ansatz (12~T|) and determine the width parameter 77 via 

t? 4 - 2v / 2^7aa 2 n 2 -1=0. (25) 

The solution of the nonlinear (1251) when substituted into (1221) yields the desired p, z [10] 
through a procedure far complicated than the analytic formulae ( 1241) . 
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The approximation scheme of MMD is quite involved but does not require the 
solution of a nonlinear variational equation like SPR. They calculate p z (n2) via [11] 

M**) = \[(JI+ VV-C 6 ) 1/3 + (i? - V^-C 6 ) 173 - C] 2 (26) 

where 77 = 4 + 6£ — ( 3 + 24nx and ( = (k — I) with 

k- 1 = ^2/^+ 6(x - 0.1)(1 - ^/7r)[l - (10x)~ 1/5 ] (27) 

where x = aa z n2, and the Heaviside step function 9(x — 0.1) = 0, for x < 0.1, and = 1 
for x > 0.1. It is realized that the function ([27)1 is not analytic in x. 



Table 2. Chemical potential [izijii) of ([20]) for different aa z ni obtained from an 
accurate numerical solution of ([H]), and from (gj) (Present, rj = 0.95), (MMD) 
P3], and (PH) (SPR) dp). 



aa z n 2 


Numerical 


MMD 


SPR 


Present 





0.5 


0.5 


0.5 


0.50263 


0.2 


1.348783 


1.31186 


1.36149 


1.34376 


1 


3.599892 


3.50165 


3.69266 


3.54355 


10 


16.45405 


16.2547 


17.0012 


16.3294 


100 


76.30080 


75.9963 


78.8855 


76.1358 



Now we calculate the chemical potential p z {n 2 ) of the three approaches and compare 
them with the precise result for p z (n2) from a numerical solution of ffT9l . The results are 
shown in table [2] for different aa z ri2 along with those from the MMD and SPR schemes. 
After a small experimentation the constant rj in (I2"4"j) was fixed at rj = 0.95 for all aa z U2- 



3. 3D-1D and 3D-2D reductions of Fermi superfluid DF equations 



We consider a Gallilei-invariant density-functional (DF) formulation for a trapped Fermi 
superfluid at BCS and unitarity limits described by [HI [7] 



-th- 
dt 



h 2 vl 

4m 



2h 2 



+ 2V{t) + 2 2/3 x |*(r 



1 4/3 



m 



V(r,t) =0, 



(2* 



with x = (3vr 2 ) 2 / 3 {/2 (£ = 1 at the BCS limit and f = 0.44 at unitarity [23]), N 
the number of fermions, and normalization f \ty\ 2 d 3 r = N/2 (|\1/| 2 is the density of 
Fermi pairs), m the mass of an atom, V(r) = mu 2 (p 2 + X 2 z 2 )/2 the harmonic trapping 
potential with frequencies u and Xu in radial (p) and axial (z) directions (f = p, z). 
The fully-paired Fermi superfluid is assumed to be composed of spin-half fermions with 
an equal number of spin-up and -down components. Employing dimensionless units 
t = icu, r = r/a p , z — z/a p , p = p/a p ,a = a/a p , ip^J N/2 = ^a 3 / 2 , a p = yjh/mu, (128!) can 
be written as 

_ Yl + p 2 + X 2 Z 2 + 2^2/3^14/31 ^ f) = Qj (2Q) 



with / |^| 2 ci 3 i 



dt 4 

1 ( I 1 2 is the density of Fermi atoms). 
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3.1. 3D-1D reduction for a cigar-shaped Fermi superfluid 



For a cigar-shaped trap, we consider the factorization ([3]). We substitute 
multiply the resultant equation by ip*(p,m) and integrate in p to get 



in $29]) and 



9 1 <9 2 - 22 . , 
«- + t— - A z - p p {n lj 



(z,t)=0, 



dt Adz 2 

where p p (rii) is the chemical potential emerging from the following 2D DF equation 



(30) 



Y7 2 

V P , 2 in 2/3 1 1 4/3 

— f + P + 2xV M 7 



yu p (n 3 



0. 



(31) 



In the small iV weak-coupling limit, the wave function </?(/?, ni) can be approximated 
by the normalized Gaussian form [21] (jSJ). With this wave form the chemical potential 
of fl3~TT) becomes 



/i p (ni) 



or 



6% 



2/3 



(32) 



4a V 5 a 4 /% 2 / 3 ' 

Chemical potential p p {ni) of f[3"2"j) is consistent with (4.7) of AS [2U]. In the large iV TF 
limit the normalization condition of the TF wave function leads to [TOl [TTj 



p p {m) 



2^J 2X 



3/5 



1.38336ni /5 x 3/5 - 



(33) 



Chemical potential (|33|) is approximately equal to the chemical potential in the large N 
limit of the corresponding model (4.10) of AS [20] which yields in present notation 

7 (6x) 3/5 n? /5 



1.42545ni /5 x 3/5 . 



(34) 



5 (57r 2 )V5 2 2/5 

Here we use the following simple interpolation formula for p p {n\) incorporating the 
limiting values (1321) and (1331 



3c^ 



4a 2 



5/3 



2X 



5ni \ 2/3 



3/5 



(35) 



where a is taken as a fixed constant for all n\. 

AS [20] solve Eq. ( I3T1) variationally and obtain for the width a 



1 

4 



12% / nia 



2/3 



(36) 

25 V 7T J K J 

Here we calculate the chemical potential p p {n\) obtained from (1351) for different a 
and from (4.7) and (4.8) of [20] and compare with the precise result for p p {ni) from a 
numerical solution of (13~T1) . Our finding is exhibited in table [3] for different ri\ together 
with those from the AS [20] scheme. After a small experimentation it was found that 
the best overall result from (1351) was found for a = 0.98/^2. 

Now to see how well the effective ID equations fl30l) and fl35l) reproduce the linear 
density 4> 2 (z) of a cigar-shaped Fermi superfluid we plot in figured] (a) the ID density 
calculated from ( 1301) and ( 135|) and the same calculated from the full 3D ( 1291 for A = 1/4 
and N = 2, 10, 100. The excellent agreement between the two sets of results for A as 
large as 1/4 demonstrates the usefulness of the present ID model equations. 
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Table 3. Chemical potential p p {ni) of (|30|) for different n\ for a Fermi superfluid at 
unitarity (£ = 0.44) obtained from an accurate numerical solution of (|31|) . from (|32]l 
and j36]) (AS) [6], and from fl35]) (Present, a = l/yg, and 0.98/\/2). 





Numerical 


AS 


Present 


Present 








a = l/y/2 


a = 0.98/\/2 





1 


1 


1 


1.00081 


0.1 


1.37401 


1.37619 


1.36784 


1.37667 


1 


2.43893 


2.46380 


2.39751 


2.41810 


10 


5.56376 


5.69170 


5.45891 


5.49147 


100 


13.70413 


14.09704 


13.5463 


13.5867 


1000 


34.29934 


35.33109 


34.1057 


34.1506 




Figure 1. (Color online) (a) The linear density 4> 2 {z) of a Fermi superfluid at unitarity 
(with £ = 0.44) vs. z (both in dimensionlcss units) calculated from the 3D DF equation 
([29]) and the ID model (JSUJ) and ([35]) for A = 1/4, a = 0.98/\/2 and N = 2, 10, and 
100. (b) The radial density (p 2 (p) of a Fermi superfluid at unitarity (with £ = 0.44) 
vs. p (both in dimensionlcss units) calculated from the 3D DF equation ([29]) and the 
2D model J37|) and flU} for A = 4, 77 = 0.92/^2 and N = 2, 10, and 100. 



3. J?. 3D-2D reduction for a disk-shaped Fermi superfluid 



In the case of a disk-shaped trap we consider the factorization (Tl4|) and substitute it in 
f[2"9"j) . multiply the resultant equation by ^*(2;,n2) and integrate in z to get 



.d 



(/?(p,t) = 0, 



0* 4 p 

where /i z (ri2) is the chemical potential emerging from the following ID DF equation 



(37) 



1 d 2 



--^ + X 2 z 2 + 2 X nT |0| 4/3 -/i,(n 2 ) 



^,n 2 ) = 0. 

It is convenient to introduce scaled variables z — z/a z ,(f> — s/a z (f), and fi z 
a z = \fljx, when ( 138|) becomes 

1 d 2 



Adz 2 



+ z 2 + 2 X a^nTW l '-Mn2) 



(z,n 2 ) = 0, 



dz(p* 



(38) 
^ 2 a^ with 

(39) 
(40) 
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The form of the chemical potential ]x z [n 2 ) of (1391) is known in the small and 
large N limits. In the small N weak-coupling limit the wave function <j>(z, n 2 ) can 
be approximated by the normalized Gaussian form [21] (T2TT) . With this wave form the 
chemical potential of (140]) becomes 

1 ^ , 2 X (n 2 al)^ [3 



2 



n 



2/3^1/3 



(41) 



Chemical potential p. z { n 2) of (j4ip is consistent with (5.7) of AS [20]. In the large N TF 
limit the normalization condition of the TF wave function leads to [TT] 

3/4 



/M«2, 



(3tt) 2 /3 



1.54947x 3/4 a 



'n 2 . 



(42) 



Chemical potential (|42p is approximately equal to the chemical potential in the large N 
limit of the corresponding model (5.10) of AS [20] which yields in present notation 

12a z 3 3 /V /4 v^ 2 



1.57383* 



3/4, 



7*.-, 



(43) 



Here we use the following interpolation formula for // z (ri2) incorporating the limiting 
values (JH]) and (1421) 



5//" 



15^V 
9vr 



/15^V 
V 9tt 



4/3 



(n 2 a 



2\2/3 



(3vr) 2 / 3 



3/4 



.(44) 



where rj is taken as a fixed constant for all a 2 z n 2 . 

AS [20] solve Eq. (1391) variationally and obtain for the width 77 

, 1 A X {a 2 z n 2 r?yi* [3 
4 5 7T X /3 V5' 



(45) 



Table 4. Chemical potential y, z [n-i) of ((40|) for different ar z ni for a Fermi superfluid 
at unitarity (£ = 0.44) obtained from an accurate numerical solution of from (|4"Tj) 
and (T431) (AS) E|, and from (Tl^ (Present, 77 = 1/V2 and 7? = 0.92/y^). 



a 2 r2. 2 


Numerical 


AS 


Present 


Present 








V = l/V2 


7] = 0.92/^ 





0.5 


0.5 


0.5 


0.50696 


0.1 


1.071781 


1.07375 


1.05892 


1.08883 


1 


2.79861 


2.82290 


2.72155 


2.78109 


10 


8.59782 


8.72108 


8.41456 


8.50454 


100 


27.09534 


27.51645 


26.8149 


26.9282 


500 


60.56659 


61.51645 


60.2317 


60.3571 



As in the 3D-1D reduction, now we calculate the chemical potential n z (n 2 ) of the 
three approaches and compare them with the precise result for £i z (n 2 ) obtained from 
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Figure 2. (Color online) (a) The rms axial size vs. time (both in dimensionless 
oscillator units) during oscillation of a cigar-shaped Fermi supcrfluid at unitarity 
(£ = 0.44) for A = 0.25, started by reducing the axial trap suddenly by a factor of 
0.9, as calculated by the full 3D equation ([29]) and the reduced ID equations (|30j) 
and (05]) for a = 0.98/-y/2- (b) The rms radial size vs. time (both in dimensionless 
oscillator units) during oscillation of a disk-shaped Fermi superfluid at unitarity for 
A = 4, started by reducing the radial trap suddenly by a factor of 0.9, as calculated by 
the full 3D equation flM]) and the reduced 2D equations ([37]) and (JM]) for rj = 0.92/V2. 



a numerical solution of (l39|) . The results are shown in table H] for different a^n 2 along 
with those from the AS scheme |20j . 

Now to see how well the effective 2D equations fl37"]) and (JHJ) reproduce the radial 
density <p 2 {p) of a disk-shaped Fermi superfluid we plot in figure [U (b) the 2D density 
calculated from (l37j) and (]44]) and the same calculated from the full 3D ([29]) for A = 4 
and N = 2, 10, 100. The excellent agreement between the two sets of results for A as 
small as 4 demonstrates the usefulness of the present 2D model equations. 

3.3. Dynamics: Breathing oscillation 

Now we subject the reduced ID and 2D DF equations to a more stringent test, e.g., 
how well these equations can reproduce non- stationary (non-equilibrium) dynamics of 
a cigar- and disk-shaped Fermi superfluid. First we consider a cigar-shaped Fermi 
superfluid with A = 0.25, which is set into breathing oscillation by reducing only the axial 
potential suddenly by a factor of 0.9. The radial trap is left unchanged. The resultant 
oscillation is studied using the full 3D DF equation ( ]29|) as well as the reduced ID DF 
equation (130]) using the chemical potential (1351) . The root mean square (rms) axial size 
as calculated from the 3D and ID equations are plotted in figure [2] (a). Next we consider 
a disk-shaped Fermi superfluid with A = 4, which is set into breathing oscillation by 
reducing only the radial potential suddenly by a factor of 0.9. The resultant oscillation 
is studied using the full 3D DF equation ( [29]) as well as the reduced 2D DF equation 
(137]) using the chemical potential (]44j) . The rms radial size as calculated from the 3D 
and 2D equations are plotted in figure [2] (b). 
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Figure 3. (Color online) (a) The rms axial size vs time (both in dimensionless 
oscillator units) during axial expansion of a cigar-shaped Fermi supcrfluid at unitarity 
(£ = 0.44) for A = 0.25, started by removing the axial trap suddenly, as calculated by 
the full 3D equation j22J and the reduced ID equations fl30|) and (35|) for a = 0.98/V2. 
(b) The rms radial size vs time (both in dimensionless oscillator units) during radial 
expansion of a disk-shaped Fermi supcrfluid at unitarity for A = 4, started by removing 
the radial trap suddenly, as calculated by the full 3D equation (|29|) and the reduced 
2D equations ([37]) and fH for -q = 0.92/^2. 

3.4- Dynamics: Free expansion 

Now we consider the problem of free expansion of a cigar- and disk-shaped Fermi 
superfluid, respectively, when the axial and radial traps are suddenly removed after 
the formation of the superfluid. First we consider a cigar-shaped Fermi superfluid with 
A = 0.25 which is allowed to expand freely in the axial direction by setting the axial trap 
suddenly to zero in the 3D equation. The radial trap is left unchanged. The resultant 
expansion is studied using the full 3D DF equation (1291) as well as the reduced ID DF 
equation ( |30l) using the chemical potential ( |35|) . The rms axial size as calculated from 
the 3D and ID equations are plotted in figure [3] (a). Next we consider a disk-shaped 
Fermi superfluid with A = 4, which is allowed to expand freely in the radial direction 
by setting the radial trap suddenly to zero in the 3D equation. The axial trap is left 
unchanged. The resultant expansion is studied using the full 3D DF equation (1291) as 
well as the reduced 2D DF equation ( 1371) using the chemical potential (1441) . The rms 
radial size as calculated from the 3D and 2D equations are plotted in figure [3] (b). The 
agreement between the dynamics as obtained from the full 3D DF equation and that 
from the reduced DF equations in Figs. [2] and [3] is quite satisfactory. 

4. Conclusion 

In conclusion, we have suggested time-dependent mean-field reduced DF equations in 
ID and 2D, respectively, for a cigar- and disk-shaped BEC and Fermi superfluid in the 
BCS and unitarity limits with simple analytic nonlinear terms. The numerical solution of 
these reduced equations reveals that they produce results for density of Fermi superfluids 
and BEC in cigar- and disk-shaped traps in excellent agreement with the solution of the 
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full 3D DF equation. We also studied non-stationary breathing oscillation of the cigar- 
and disk-shaped Fermi superffuid initiated by a sudden change of axial and radial traps, 
respectively Finally, we applied the reduced equations to the study of free expansion 
of a cigar- and disk-shaped Fermi superfluid initiated by a sudden removal of the axial 
and radial trap, respectively. The reduced equations produced equally good results in 
both these studies when compared with the solution of the full 3D equations. 
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